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The second J 77 = 2 + state of 12 C, predicted over fifty years ago as an excitation of the Hoyle 
state, has been unambiguously identified using the 12 C(7, ao) 8 Be reaction. The alpha particles 
produced by the photodisintegration of 12 C were detected using an Optical Time Projection Cham- 
ber (O-TPC). Data were collected at beam energies between 9.1 and 10.7 MeV using the intense 
nearly mono-energetic gamma-ray beams at the HI7S facility. The measured angular distributions 
determine the cross section and the E1-E2 relative phases as a function of energy leading to an 
unambiguous identification of the second 2 + state in 12 C at 10.03(11) MeV, with a total width of 
800(130) keV and a ground state gamma-decay width of 60(10) meV; B(E2: 2+ -> 0+) = 0.73(13) 
e 2 fm 4 [or 0.45(8) W.u.]. The Hoyle state and its rotational 2 + state that are more extended than 
the ground state of 12 C presents a challenge and constraints for models attempting to reveal the 
nature of three alpha particle states in 12 C. Specifically it challenges the ab-initio Lattice Effective 
Field Theory (L-EFT) calculations that predict similar r.m.s. radii for the ground state and the 
Hoyle state. 



The second J w = 0+ state at 7.654 MeV in 12 C, first 
predicted by Hoyle [lj in 1953 and thus called the Hoyle 
state, plays a central role in nuclear physics. It is a well 
known fundamental testing ground of models of the clus- 
tering phenomena in light nuclei which is highlighted by 
recent developments of ab-initio theoretical calculations 
that are able to calculate light nuclei such as 12 C. The 
Hoyle state plays a central role in stellar helium burn- 
ing by enhancing the production of 12 C in the universe 
allowing for life as we know it. It is the first and quite 
possibly still the best example of an application of the 
anthropic principle in physics. Early on after the discov- 
ery of the Hoyle state it was suggested by Morinaga [|[ 
that we can learn more about the structure of the Hoyle 
state by studying the rotational band built on top of it 
which led to a fifty year long search for the second 2 + 
state in 12 C 0. 

Recently the existence of the second 2 + state in 12 C has 
been the subject of much debate. It was observed at ap- 
proximately 9.8 MeV in measurements of the 12 C(ev,cv') 
and the 12 C(p,p') inelastic scattering reactions 0-0]) bvrt 
it was not observed below 10 MeV in either the beta- 
decays of 12 N and 12 B [8[ or in a recent measurement 
of 3 He induced reactions on 10 > n B at 4.9 and 8.5 MeV 
0. In contrast, analysis of the beta-decay data @ sug- 
gest a 2 + state at 11.1 MeV which was not observed in 
the 12 C(p,p') data @ or the recent measurement of the 
n B( 3 He,d) reaction at 44 MeV [l|. 

Previous measurements @, H, 0, @| are dominated by 
the broad (r « 3.0 MeV) 0+ state at 10.3 MeV, as 
well as the narrow 3" state at 9.641 MeV [Hi]. Indeed 



a 2 + state below 10.0 MeV was observed in the inelas- 
tic scattering data [J, Q only after separating from a 
large background contribution from the third + at 10.3 
MeV. Gamma-ray beams as used in this study [TTf] can- 
not populate + states and will populate the 3~ state 
with very small probability making them an excellent 
probe to use in the search for 2+ state in 12 C. Since 7- 
ray beams can also induce El transitions the nearby 1~ 
state at 10.84 MeV is expected to contribute. In the cur- 
rent study we used the Optical Time Projection Chamber 
(O-TPC) detector discussed in [l^] to detect the outgoing 
particles with nearly 100% efficiency. Thus our gamma- 
ray beam plus an O-TPC detector system is a powerful 
(beam-target-detector) combination in the search for 2 + 
states in 12 C. In this paper we present background free 
data with an unambiguous identification of the second 
2+ state at 10.03(11) MeV in 12 C. 

Ever since Brink suggested that the Hoyle state is a 
very extended object with the structure of three alpha- 
particles arranged in a linear chain [131 ] many theoretical 
models have been developed to describe the Hoyle state 
and the structure of 12 C. One of the issues of great cur- 
rent interest is the r.m.s. radius of the Hoyle state and in 
particular whether the Hoyle state is an extended object 
with an rms radius considerably larger than the ground 
state of 12 C. 

A number of models have been proposed to describe 
the structure of 12 C including an algebraic U(7) model 
with a D3/j symmetry [l4l ]. a microscopic Fermionic 
Molecular Dynamic (FMD) model [l5| together with a 
"BEC-like" cluster model [16[ which predicts wave func- 
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tions that are quite similar to those previously pre- 
dicted by the RGM cluster model [l7|, an ab-initio No- 
Core Shell Model (NCSM) [3, and a No-Core Simplec- 
tic Model (NCSpM) and ab-initio Lattice Effective 
Field Theory calculations (L-EFT) [20|, Sj. The NCSM 
calculations that currently extend up to 10 fiuj do not ad- 
equately predict the location of the Hoyle state [l8| and 
in the symmetry inspired schematic NCSpM calculations 
fl9| a model space of up to 20 Twj is needed to predict the 
low lying Hoyle state close to the measured energy. The 
many fwj model space required in these shell model calcu- 
lations is suggestive of a very extended alpha-clustering 
configuration. Such an extended state arises naturally in 
cluster models [14T417I ] and it represents a major challenge 
to ab-initio calculations [lg, 120, |2l| . 

Current models differ on the shape of the Hoyle state. 
In the U(7) model [l| and the FMD model [H the Hoyle 
state is predicted to be an oblate equi-lateral triangular 
three alpha-particle configuration. Both U(7) and FMD 
models predict a rotational band built on the Hoyle state 
but the FMD model predicts the second 2 + in 12 C to have 
a B(E2) to the third + state which is twice as large as 
that leading to the Hoyle state, and thus not to be a 
member of the Hoyle rotational band. In the L-EFT 
calculations [2l[ the Hoyle state is primarily of the bent- 
arm chain (or obtuse triangular) shape. Both NCSpM 
calculations [l9[ and L-EFT calculations !2l| predict the 
Hoyle state to have a deformed prolate shape with a ro- 
tational band built on it. The "BEC-Like" cluster model 
predicts the Hoyle state to be spherically symmetric. 
In addition while the FMD model predicts an rms radius 
of the Hoyle state (3.38 fm) that is larger than the 
ground state (2.34 fm) [15|, the ab-initio L-EFT calcula- 
tions [2l| predicts an rms radius (2.4 fm) equal (within 
the predicted error bar) to the rms radius of the ground 
state of 12 C. The NCSpM calculations [H| predict an rms 
radius (2.93 fm) that is 25% larger than the ground state 
of 12 C. 

The three alpha-particle structure of 12 C naturally 
leads to models that utilize triangular geometry [3, [T^, 
US [HI- Such triangular systems are ubiquitous in physics 
including the Xj molecular system [22| and the three 
quark system [23, H3] , and their spectra resemble the one 
predicted by the oblate spinning top with a D^h symme- 
try [H, [23)]. Unlike molecules, in nuclei the energy scale 
of rotations and vibrations are similar, leading to large 
mixing of rotational and vibrational states which leads to 
deviation from the prediction of a rigid rotor jli| . Still 
the phenomenological schematic U(7) model preserves 
the rotation-vibration structure [L4j and it serves as a 
useful guiding tool for discussing the essential degrees of 
freedom of the three alpha-particle system. 

The current measurement of the 12 C(7, a) 8 Be reac- 
tion was performed at the HI7S facility that produces 
an intense, nearly monoenergetic gamma-ray beam by 
Compton backscattering photons of free-electron laser 
[lll |. Beams of circularly polarized gamma-rays with 
energies between 9.1 and 10.7 MeV were used with en- 




FIG. 1: (Color Online) A typical image recorded by the 
CCD camera of three alpha-partciles from the reaction 
12 C( 7 ,a ) 8 Be a + a). 



ergy spreads of 300 - 350 keV and on-target intensities 
of w 2 x 10 8 7/sec. The beam intensity was measured 
by detecting neutrons from the d(7,n)p reaction using 
an in-beam D2O target, cross-calibrated against a large 
Nal(Tl) detector. The energy profile of the beam was 
measured using a large HPGe detector, and the spectra 
were unfolded using a Monte Carlo technique (2|| [26| . 
The alignment of the detector with respect to the beam 
was achieved using a gamma camera and lead absorbers 
placed in the front and back of the detector as discussed 
in 0. 

An O-TPC operating at 100 torr with the gas mix- 
ture of CO 2 (80%) + N 2 (20%) [13 was used to detect 
the outgoing alpha-particles from the 12 C(7, a) 8 Be re- 
action. The events recorded in the O-TPC include pro- 
tons from the 14 N(7,p) reaction, alpha particles from the 
16 ' 18 0(7,a) and the 12 C(7, a) 8 Be reactions, and cosmic 
rays. Nearly all (98%) of the 12 C dissociation events 
were observed to proceed via the 12 C(7,ao) 8 Be reaction 
leading to the ground state of 8 Be, and the subsequent 
immediate decay to two nearly co-linear alpha-particles 
as shown in Fig. 1. The 12 C(7, a) 8 Be events were easily 
separated from the other events listed above except for 
the 16 0(7, a) events, using the energy deposited in the 
detector (with a measured energy resolution of FWHM 
« 100 keV [Hj])- The recorded energy and track from the 
12 C and 16 O dissociation events are very similar but the 
events can be distinguished using the line-shape analysis 
of the time projection signals recorded by the photomul- 
tiplier tubes jl2l ]. 

Each measured (PMT) time projection signal was fit 
using the calculated line shapes of the 12 C(7,a) 8 Be and 
16 0(7,a) 12 C events [Hj]. A good fit of all 12 C and 16 
events was obtained. The goodness-of-fit parameters, Xc 
and Xo °f the predicted line shapes of 12 C(7, a) 8 Be and 
16 0(7, a) 12 C respectively, were used to classify the events 
as shown in Fig. 2 for the beam energy of 9.8 MeV. All 
events to the left of the dotted (red) line were identified 
as 12 C(7, a) 8 Be events. The efficiency of the cut and 
leakage of 16 0(7, a) 12 C events were estimated by fitting 
the distribution to the sum (shown by solid black line) of 
two log- normal functions (shown by dashed blue lines). 
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FIG. 2: (Color Online) Event identification function (E 7 = 
9.8 MeV) derived from the goodness-of-fit parameters (xc 
and Xo f° r 12 C an d 16 dissociation events, respectively) as 
discussed in the text. 



The cut was placed such that fewer than 0.5% of the 
events to the left of the cut were 16 0(7, a) 12 C events. 

Complete angular distributions of 12 C (7, a ) 8 Be 
events were measured at seven energies between 9.1 and 
10.7 MeV. The events recorded in the O-TPC are trans- 
formed to the (6, <fi) coordinate system [12] with an accu- 
racy in 9 varying between 2.5° and 6.0°, depending on the 
out-of-plane angle of the track. The angular distributions 
were fit in terms of El + E2 amplitudes and their relative 
phase 4>i2 as discussed in section 4.1 of 27] . Since angular 
information was available for each 12 C(7, ao) 8 Be event 
individually, unbinned maximum likelihood fits were used 
to avoid losing information through binning. Angular 
distributions measured at gamma-ray beam energies of 
of 9.6 and 10.7 MeV are shown in Fig. 3 along with the 
fit yielding the cross section ratio of EiMl = 0.97tp 02 
and 0.711qq 5 , as well as phase angles </>i 2 = 80 ± 6° and 
132 ± 5°, respectively. The angular distributions were 
dominated by the E2 component at all but the highest 
beam energy (at 10.7 MeV) where a non-negligible con- 
tribution of the 1~ state at 10.84 MeV leads to a very 
asymmetric angular distribution, as shown in Fig. 3. 

The total El and E2 cross sections and the relative El- 
E2 phase ((^12) extracted using the angular distribution 
data (as shown in Fig. 3) are shown in Fig. 4 as a func- 
tion of energy with error bars that include both statistical 
and systematic uncertainties. The systematic uncertain- 
ties associated with each measured cross section are dom- 
inated by a 5% uncertainty in the gamma-ray beam in- 
tensity. In Fig. 4(a) we show the El and E2 cross section 
components measured at these energies together with fits 
to Brcit-Wigner resonances with energy-dependent level 
shifts and widths 28], convoluted with the measured 
gamma-ray beam energy distribution. Coulomb wave 
functions were calculated using the continued-fraction ex- 
pansion technique [29j j with tq — 1.4 fm. The fit to the 
El cross section data uses the previously determined en- 
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FIG. 3: (Color online) Angular distribution for 12 C(7, a ) 8 Be 
events measured at a beam energy of 9.6 and 10.7 MeV. The 
solid curve is the fit that included El and E2 amplitudes as 
discussed in the text. The error bars are statistical only. 
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FIG. 4: (Color online) (a) The measured El and E2 cross 
sections of the 12 C(7, ao) 8 Be reaction, (b) The measured El- 
E2 relative phase angle (^12) together with the phase angle 
calculated from a two-resonance model. 



ergy and width of the 1" state at 10.84 MeV in 12 C pj. 
with the strength adjusted to fit the data. The fit to 
the E2 data includes three free parameters: the partial 
widths (r a , r 7 ) and the resonance energy. 

The E2 cross section data allow us to identify a 2 + 
resonance at 10.03(11) MeV with a total (alpha-particle) 
width of 800(130) keV (that exhausts 65(9)% of the 
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Wigner limit) in agreement with Ref [7|. The measured 
gamma-decay width to the ground state is 60(10) meV 
leading to a reduced quadrupole electromagnetic transi- 
tion of B(E2 : 2+^0^) = 0.73(13) e 2 fm 4 [or 0.45(8) 
W.u.]. This measured B(E2) is not too different from the 
prediction of the FMD model (0.46 e 2 fm 4 ) [H], but some- 
what smaller than predicted in the L-EFT calculations 
[2(1) e 2 fm 4 ] [2l|. Note that the slight difference between 
the maximum of the calculated cross section (at 9.8 MeV) 
and the resonance energy (at 10.03 MeV) is due to the 
energy-dependent widths used in the fit which push the 
maximum yield toward lower energies. The highest en- 
ergy data point at 10.7 MeV seems inconsistent with this 
single resonance. In order to estimate the error in the 
measured resonance energy, we also analyzed our data 
including another 2 + that was previously suggested at 
11.2 MeV Q leading to a total error in the resonance 
energy of 110 keV. The current results do not allow us to 
place any constraints on 2 + states at energies above 11 
MeV. 

The E1-E2 phase differences ((^12) extracted from our 
measured angular distributions are shown in Fig. 4(b). 
The measured phase angle is compared to the predicted 
phases [27j : 

<fii2 = 62 — Si + arctan(r]/2) 

where the nuclear phase shifts 8e are given by the reso- 
nance phase shift minus the hard sphere phase shift [28j 
and 77 is the Sommerfcld parameter. The I = 1 reso- 
nance phase shifts (8\) were calculated using parameters 
of the known l - state at 10.84 MeV [3(| and the dip is 
due to the hard sphere contribution. The I = 2 phase 
shifts {82 ) were calculated using the resonance energy 
and width determined from the fit to the E2 cross sec- 
tion data. The calculated </>i2 curve was averaged over 
the measured cross section and gamma-ray beam energy 
distribution. The good agreement between the measured 
and calculated phases unambiguously establishes the ex- 
istence, the energy and the width of the 2 + resonance 
reported here and indicates that there is little or no con- 



tribution from other (background) amplitudes. 

The measured second 2+ state at 10.03(11) MeV re- 
ported in this work lies 2.38(11) MeV above the + Hoyle 
state which is approximately half the excitation energy 
of the first 2 + state of 12 C. Since the U(7) model predicts 
the ground state rotational band and the Hoyle bands to 
arise from the same geometrical shape we conclude that 
in this model the radius parameter of the Hoyle state is 
approximately \/2 larger than the r.m.s. radius of the 
ground state. This conclusion is consistent with (but 
slightly larger than) the r.m.s. radius of the Hoyle state 
[2.89(4) fm] determined from inelastic light- ion scattering 
experiments of [3l[. It should be noted that the present 
result and all previous determinations of the rms radius 
of the Hoyle state using 12 C(e,e') data [l5[ and 12 C(x,x') 
data |31j are model dependent. Our inability to measure 
elastic scattering off the (very short lived) Hoyle state 
makes a direct measurement of the r.m.s. radius of the 
Hoyle state unlikely. 

In conclusion we have used an Optical Time Pro- 
jection Chamber (O-TPC) with intense nearly mono- 
energetic gamma-ray beams to unambiguously identify 
the long sought after second 2 + state in 12 C having an 
extended three alpha-particle configuration. We provide 
the resonance parameters including the spin and parity 
( J n = 2+), energy, total (alpha-particle) width, and the 
B(E2) value for the decay to the ground state which must 
be reproduced by theoretical models which properly de- 
scribe the low lying three alpha-particle structure of 12 C. 
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